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ABSTRACT 

Relaxation  of  the  density  operator  of  an  atom  near  a  four-wave  mixing 
phase  conjugator  (PC)  is  studied.  The  relaxation  operator  F  can  be  expressed 
in  the  Cartesian  components  of  the  atomic  dipole  operator  and  correlation 
functions  of  th  /acuum  electromagnetic  field.  We  evaluate  these  correlation 
functions  in  terms  of  the  Fresnel  coefficients  for  reflection  and  transmission 
of  a  polarized  monochromatic  plane  wave.  Our  expression  for  F  includes  both 
linear  and  nonlinear  interactions  in  the  medium,  and  it  reduces  to  well-known 
results  in  the  limits  of  a  pure  dielectric,  a  mirror,  and  empty  space.  The 
example  of  a  model  two-state  atom  in  combination  with  a  model  PC  is  worked  out 
in  detail.  From  the  rate  equations  for  the  populations  of  the  levels  we  infer 
that  the  relaxation  acquires  a  contribution  from  both  ordinary  spontaneous 
decay  and  stimulated  transitions  which  are  induced  by  the  PC.  It  is  shown 
that  an  atom  in  its  ground  state  has  a  finite  probability  of  being  excited, 
and  that  the  excited  state  has  a  finite  population  in  the  long-time  limit.  We 
also  work  out  the  case  of  a  degenerate  two-level  atom  in  combination  with  a 
realistic  model  of  a  transparent  PC,  and  it  is  shown  that  our  present  results 


2 


for  the  transition  rates  are  consistent  with  earlier  calculations  of  the 
fluorescence  yield. 


PACS;  42.65.H,  32.80,  52.50.D,  42. 50. Kb 
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I.  INTRODUCTION 

An  excited  atom  can  decay  spontaneously  to  a  lower  state,  in  combination 
with  the  emission  of  a  fluorescent  photon,  due  to  the  coupling  of  its  dipole 
moment  p  with  the  electromagnetic  field  of  the  vacuum.  For  a  model  two- state 
atom  with  excited  state  j  e>  and  ground  state  lg>,  and  level  separation 
this  amounts  to  an  effective  lifetime  1/A  of  the  excited  state,  where 

3 

w  2 

A  -  ^  |<e|M|g>|  ,  (1.1) 

3irc  ]iic 
o 

equals  the  Einstein  coefficient  for  spontaneous  decay.  This  relaxation 
constant  A  pertains  to  an  atom  in  empty  space.  When  an  atom  is  close  to  the 
surface  of  a  dielectric  or  a  mirror,  with  the  atom-surface  normal  distance  of 
the  order  of  an  optical  wavelength  (which  is  large  compared  to  the  atomic 
dimensions),  then  the  expression  for  the  lifetime  of  an  excited  state  changes 
dramatically  [1-12].  Such  an  alteration  of  lifetimes  due  to  the  presence  of 
media  was  first  demonstrated  by  Drexhage  for  optical  transitions  in  dye 
molecules  in  thin  films  [13].  Much  later,  lifetime  measurements  were 
performed  on  a  single  Rydberg  atom  in  a  microwave  cavity  [14,15],  and  for  an 
ordinary  atom  near  a  mirror  or  in  between  parallel  mirrors  Ho,!?].  It 
appeared,  in  agreement  with  theory,  that  both  an  enhancement  a  .d  an  inhibition 
of  the  spontaneous  decay  rate  could  be  obtained. 

In  order  to  understand  this  phenomenon  one  merely  has  to  realize  that  a 
photon,  which  is  emitted  by  the  atom  into  the  direction  of  the  surface,  has  a 
probability  of  being  reflected  and  a  probability  of  being  refracted,  just  as 
in  classical  optics.  A  reflected  wave  then  interferes  with  photons  which  are 
emitted  in  a  direction  away  from  the  medium,  and  this  can  give  rise  to  either 
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a  reduction  or  an  increase  of  the  fluorescent  intensity.  Since  energy 
conservation  requires  that  the  energy  which  is  emitted  by  the  atom  as 
radiation  must  be  provided  by  that  atom,  we  conclude  that  a  destructive 
(constructive)  interference  of  radiation  must  be  accompanied  by  an  increase 
(decrease)  of  the  lifetime  of  the  atomic  state.  From  a  different  point  of 
view,  we  can  say  that  the  electromagnetic  vacuum  near  a  medium  is  different 
from  the  vacuum  of  empty  space.  Without  the  medium,  the  electric- field 
operator  in  the  Heisenberg  picture  is  a  superposition  of  travelling  plane 
waves,  without  any  preference  for  a  specific  orientation  in  space.  An  excited 
atom  can  emit  a  fluorescent  photon  in  any  of  the  available  plane-wave  modes, 
and  these  photons  travel  away  from  the  atom  (their  superposition  being,  of 
course,  spherical  dipole  radiation).  Near  an  interface,  however,  this  mode 
structure  is  different,  and  therefore  the  coupling  of  the  atomic  dipole  fi  to 
this  vacuum  field  will  be  different.  Since  this  interaction  is  responsible 
for  the  spontaneous  decay  and  the  emission  of  fluorescence,  the  atom  will 
behave  differently  near  a  medium.  For  instance,  the  atom  cannot  excite  a 
plane-wave  mode  of  the  radiation  field  with  a  wave  vector  in  the  direction  of 
the  medium,  without  also  partially  exciting  the  mode  with  its  wave  vector  in 
the  specular  direction. 

We  consider  an  atom,  with  dipole  moment  fi,  which  is  positioned  at  r  - 
he^,  h  >  0,  on  the  z-axis,  and  the  plane  z  -  0  separates  the  vacuum  in  z  >  0 
from  a  nonlinear  medium  in  z  <  0  (not  necessarily  half -  infinite) .  Two  strong 
counterpropagating  laser  beams  with  frequency  w  >  0  pump  the  medium  (they 
travel  parallel  to  the  z  -  0  plane),  such  that  a  four-wave  mixing  process 
between  the  pumps  and  an  external  field,  which  is  incident  from  the  region 
z  >  0,  is  responsible  for  the  reflection  of  radiation  back  into  the  region  z  > 
0.  This  'device'  is  called  a  phase  conjugator  (PC).  We  shall  not  assume  any 
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details  about  the  four-wave  mixing,  nor  about  the  linear  properties  of  the 
medium,  which  keeps  the  theory  very  universal.  It  can  be  shown  in  general 
[18]  that  when  a  plane  wave  with  wave  vector  k  and  frequency  w  >  0  is  incident 
upon  a  PC,  then  there  are  two  plane  waves  reflected  back  into  the  z  >  0 
region.  First  there  is  of  course  the  specular  wave  with  wave  vector  k^  and 
frequency  cj  >  0 ,  due  to  a  dielectric  constant  which  is  unequal  to  unity. 
This  radiation  (when  it  comes  from  the  dipole)  is  responsible  for  a  change  in 
lifetime  of  excited  states,  as  described  above.  In  addition,  there  is  a 
phase -conjugated  wave  with  wave  vector  k^^  =  k,  and  this  wave  has  the  negative 
frequency 

u)'  —  u>  -  2u)  .  (1-2) 

Hence  the  phase  -  conj  ugated  (pc)  wave  and  the  incident  wave  are  almost 
counterpropagating,  in  contrast  to  the  propagation  direction  of  the  specular 
wave  (which  we  shall  call  the  r-wave). 

For  pure  plane  waves  it  does  not  seem  to  make  much  of  a  difference  in 
which  direction  the  produced  wave  emanates  from  the  medium:  it  is  again  a 
plane  wave.  Dipole  radiation  is  a  spherically-diverging  wave,  and  for 
reflection  at  a  linear  medium  (as  r-waves)  the  total  reflected  field  is  again 
a  diverging  spherical  wave.  The  phase -conjugated  field  (pc-waves),  however, 
is  a  converging  wave  which  is  focused  almost  exactly  on  the  radiating  atom, 
because  every  Fourier  plane-wave  component  retraces  its  path  after  reflection. 
Most  r-wave  photons  are  scattered  away  by  the  medium,  but  every  pc-wave  photon 
returns  to  the  atom  and  can  again  interact  with  it.  Such  a  process  can 
therefore  be  anticipated  to  give  rise  to  much  stronger  effects  on  atomic 
lifetimes  and  on  the  time  evolution  of  the  atomic  density  operator,  as 
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compared  to  the  influence  of  r-wave  photons.  As  it  turns  out,  however,  the 
situation  is  more  subtle,  and  simple  interpretations  as  given  above  for  linear 
reflection  have  to  be  modified  carefully  in  order  to  picture  a  consistent 
interpretation  of  the  behavior  of  an  atom  near  a  PC. 

II.  RELAXATION 

The  interaction  of  the  atomic  dipole  ^  with  the  electromagnetic  field  of 

the  vacuum  gives  rise  to  both  spontaneous  decay  and  the  emission  of 

fluorescence,  and  energy  conservation  obviously  requires  a  relation  between 

both  processes.  Elsewhere  [19]  we  have  derived  an  expression  for  the 

fluorescence  radiation  field,  as  it  is  emitted  by  the  atom  near  the  PC,  and 

the  result  was  expressed  in  terms  of  fi(t) ,  where  the  t-dependence  signifies 

the  Heisenberg  picture.  In  evaluating  quantum  expectation  values  of  certain 

field  properties  (like  the  intensity  of  emission,  or  the  spectral  distribution 

of  the  fluorescence),  we  can  always  transform  to  the  Schrbdinger  picture  in 

which  the  time  evolution  is  governed  entirely  by  the  time  dependence  of  the 

atomic  density  opertaor  p^(t) .  Spontaneous  decay  amounts  to  a  relaxation  of 

p  (t),  and  in  this  way  is  the  temporal  evolution  of  the  radiation  field 
a 

determined  by  the  time  dependence  of  the  atomic  density  operator.  In  this 
paper  we  study  the  equation  of  motion  for  p^(t),  and  a  combination  with  the 
results  from  our  previous  paper  then  completely  determines,  in  principle,  the 
temporal  and  spectral  properties  of  phase-conjugated  fluorescence. 

Since  the  atom  and  the  radiation  field  interact,  we  have  to  describe 
their  joint  quantum  state  by  a  single  density  operator  p(t).  However,  we  are 
only  interested  in  the  state  of  the  atom,  irrespective  of  the  state  of  the 
field.  This  atomic  state  is  defined  by 


* 
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p  (t)  -  Tr  p(t)  .  (2.1) 

a  L 

where  the  trace  runs  over  the  states  of  the  field  only.  In  general,  the  time 
evolution  of  p(t)  is  given  by 

ix  ^  -  [H.p]  ,  (2.2) 

and  the  Hamiltonian  has  the  form 

H  -  H  +  H  +  H  ,  (2.3) 

a  r  ar 

in  obvious  notation.  Then  it  is  a  standard  procedure  [20-22]  to  derive  an 
equation  of  motion  for  p  (t)  from  Eq .  (2.2),  under  the  condition  that  the 
radiation  field  can  be  considered  as  a  thermal  reservoir,  which  is  no 
approximation  at  all  for  the  present  problem.  Also,  we  make  the  Markov 
approximation,  which  is  perfectly  justified  for  the  interaction  between  a 
dipole  and  the  vacuum  radiation  field.  We  shall  adopt  the  compact  Liouville 
notation,  since  this  appears  to  be  the  most  convenient  framework  to  describe 
relaxation  phenomena  [23].  With  the  definition  of  the  Liouvillians 

-  X  ^  [Hj^,c7]  ,  i  -  a,r,ar  ,  (2.4) 

where  a  is  an  arbitrary  operator  in  Hilbert  space,  the  equation  of  motion 
reads 


L  +  L  )p 
r  ar 


(2.5) 
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From  Eq. 


(2.5)  one  then  derives  the  equation  of  motion  for  p  (t) , 

a 


'  iT-  •  (2.6) 


where  the  relaxation  operator  F  is  found  to  be 


Fa 

a 


-i(L  +L  )t 
a  r 
e 


L 

ar 


iL  r 

e  ^  (a  p  ) 
a  r 


(2.7) 


This  relation  defines  the  action  of  F  on  an  arbitrary  atomic  operator  ,  and 
it  involves  the  thermal - equi 1 ibrium  density  operator  p^  of  the  radiation 
field.  For  the  present  situation  we  have 


Pp  -  |0X0|  .  (2.8) 

which  Indicates  that  the  electromagnetic  field  is  in  the  vacuum  state.  In  Eq. 

(2.6) ,  the  Liouvillian  L  represents  the  free  evolution  of  the  atom,  as  it 
would  be  without  the  coupling  to  the  radiation  field,  and  the  relaxation 
operator  F  accounts  for  spontaneous  decay.  Its  matrix  elements  are  related  to 
the  lifetimes  of  atomic  states,  and  they  include  the  time  evolution  of  atomic 
coherences.  The  purpose  of  this  paper  is  to  evaluate  F,  as  it  is  given  by  Eq. 

(2.7) ,  for  an  atom  near  a  PC,  and  to  interpret  the  physical  behavior  of  the 
atom  from  the  solution  of  Eq.  (2.6). 

III.  DIPOLE  COUPLING 

In  order  to  evaluate  expression  (2.7)  for  F,  we  need  to  specify  the 

interaction  Hamiltonian  H  We  write  E  (r,t)  for  the  electric-field  operator 

ar  -V  - 
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A 


of  the  vacuum,  and  we  take  the  Schrodinger  and  Heisenberg  representation  to 
coincide  at  t  -  0.  With  h  the  position  of  the  atom,  the  interaction 
Hamiltonian  in  the  dipole  approximation  then  reads 


H  -  (h,0) 

ar  c  _v  - 


(3.1) 


The  explicit  form  of  E^(r,t)  is  irrelevant  for  the  remainder  cf  this  section. 
When  we  substitute  expression  (3.1)  into  Eq .  (2.7)  and  work  out  the 

commutators  which  define  the  various  Liouvillians ,  then  it  appears  that  F  can 
be  represented  by  the  simpler  form 


i 


(3.2) 


where  the  summation  runs  over  the  Cartesian  coordinates  i  -  x,y,z.  Here  we 
have  dropped  the  subscript  'a'  on  the  atomic  operator  o,  as  we  shall  do  in  the 
remainder  of  the  paper.  Also,  we  shall  write  p  for  the  density  operator, 
instead  of  p  The  Hilbert  -  space  operator  Q.  (not  a  Liouvillian)  is  defined 

S.  1 

as 


^i 


iL 


f .  .  (Oit. 

U  J 


(3.3) 


in  terms  of  the  nine  functions  f^j(r).  These  quantities  are  the  vacuum 
correlation  functions,  which  are  given  by 


(3.4) 
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and  where  v'e  have  used  Eq.  (2.8)  for  Notice  that  the  right-hand  side  of 

Eq.  (3.2)  does  not  contain  field  Liouvillians  anymore,  and  that  the  operators 


Q.  only  involve  the  Liouvi Ilian  L  of  the  free  evolution  and  the  Cartesian 
components  of  the  dipole  operator.  All  necessary  information  about  the 
vacuum  radiation  field  is  contained  in  the  nine  correlation  functions  f.,(r). 


IV.  PLANE-WAVE  MODES 

Before  we  write  down  the  expression  for  the  vacuum  field  ^  *'r,t),  we 
briefly  summarize  the  fundamental  plane-wave  mode  solution  for  a  PC,  mainly  to 
set  up  the  notation.  The  nonlinear  medium  is  assumed  to  occupy  the  region 
0  >  z  >  -A,  with  A  >  0  the  layer  thickness.  Then  we  consider  a  monochromatic 
polarized  plane-wave  as  the  incident  field  from  the  region  z  >  0,  with 
positive- frequency  part 


E.  (r,t) 

-me  - 


(  +  ) 


-  E,  e,  e 
Ko-ka 


i(k»r-wt) 


(4.1) 


Here,  w  -  ck  >  0  and  k  is  assumed  to  be  real  (travelling  wave).  For  later 
purposes  we  allow  the  amplitude  factor  E^  to  be  a  quantum  operator.  For  the 
unit  polarization  vector  e^^^  we  take  either  a  -  s  or  a  -  p ,  corresponding  to  a 
surf  ace  -  polarized  and  a  plane -polarized  wave,  respectivley .  Then  it  can  be 
shown  [18,24]  that  the  reflected  waves,  the  transmitted  waves,  and  the  waves 
inside  the  medium  are  also  plane  waves,  and  they  have  the  same  polarization  a 
as  the  incident  wave.  Also,  when  we  write  for  the  wave  vector  of  the  incident 
field 


k  “  ki|  +  k  e 
-li  z-z 


(4.2) 


where  the  subscript  ||  refers  to  the  plane  z  -  0,  then  every  other  wave  must 
have  a  wave  vector  with  the  same  parallel  component  k||  ,  due  to  the  boundary 
conditions  at  z  -  0  and  z  -  -A.  Therefore,  every  wave  vector  can  be  written 
as 


k. 

-1 


+  k.  e 
1  ,  z-z 


(4.3) 


and  the  value  of  k^  ^  is  determined  (up  to  its  sign)  by  the  dispersion 

relation  for  either  the  medium  or  the  vacuum.  For  a  wave  in  vacuum  with  wave 

vector  k.  we  take  the  unit  polarization  vectors  as 
-1 


X  e 
~z 


(4.4) 


X 


(4.5) 


The  waves  in  the  medium  are  in  general  not  transverse,  and  the  expressions  for 
their  polarization  vectors  are  more  complicated  [24].  For  the  present 
probl^;m,  however,  we  need  only  the  waves  in  vacuum,  so  that  we  shall  not  write 
down  the  form  of  the  fields  in  the  medium. 

For  given  values  of  k  -  w/c  and  k||  (or  the  angle  of  incidence),  the  z- 
component  of  the  wave  vector  of  the  incident  field  is 


(4.6) 


As  mentioned  in  the  Introduction,  there  are  in  general  two  waves  reflected 


back  into  the  region  z  >  0,  and  they  have  wave  vectors  k^  and  k 


From  the 


L2 


dispersion  relation,  and  the  fact  that  both  waves  must  travel  into  the 
positive  z-direction,  we  find  the  z-components  of  these  wave  vectors  to  be 


(4.8) 

where  we  introduced  the  dimensionless  detuning  parameter 


p  —  -W'  /iji 


(4.9) 


in  terms  of  w'  from  Eq.  (1,2).  Then  the  positive-frequency  part  of  the  total 
field  in  z  >  0  assumes  the  form 


i(k»r-wt)  i(k  ^r-wt) 

-  -r 


.  -i(k  ‘r+pwr) 

+  E^  P*  e  e  ~ 

k<7  k<7-k  a 

-  -  -pc 


(4.10) 


in  terms  of  the  Fresnel  reflection  coefficients  R,  and  P,  for  the  r-wave  and 
the  pc-wave,  respectively.  These  Fresnel  coefficients  contain  detailed 
information  about  the  linear  interaction,  the  four-wave  mixing  process,  the 
layer  thickness,  etc.  They  can  be  found  by  solving  Maxwell's  equations  for 
the  given  incident  wave,  although  the  result  depends  in  a  very  complicated  way 
on  the  various  parameters  [18,24],  It  is  important  to  notice  the  dagger  in 
the  last  term  on  the  right-hand  side  of  Eq .  (4.10).  For  an  incident  field 
with  frequency  o)  >  0,  the  pc-wave  is  emitted  at  the  negative  frequency  ,  and 
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which  follows  from  the  Hermiticity  of  E,  we  then  obtain  a  dagger  on  (and  a 

star  on  in  Eq .  (4.10).  For  a  classical  field  this  is  a  complex 

conjugation,  which  is  where  the  term  'phase  conjugation'  originates  from. 

Notice  that  also  the  spatial  part  exp(-ik  ‘r)  =  exp(-ik*r)  is  the  complex 

—pc  —  —  — 

conjugate  of  the  spatial  part  of  the  incident  field,  apart  from  a  small  shift 
in  wave  number  (due  to  p  1)  .  In  quantum  mechanics  this  complex  conjugation 
translates  into  a  Hermitian  conjugation,  as  is  well  known  [25-27]. 

The  incident  wave  travels  through  the  layer  and  exits  the  medium  at  the 
interface  z  -  -A.  This  transmitted  wave  (t-wave)  has  wave  vector  k^,  and 
obviously  we  must  have 

k^  -  k  ,  ('!^.12) 


whereas  the  frequency  of  the  wave  equals  w.  When  this  incident  wave  hits  the 
boundary  z  -  -A,  then  a  part  of  it  will  be  reflected  back  into  the  medium,  and 
approximately  into  the  specular  direction,  just  as  in  linear  optics.  Inside 
the  medium  this  positive  -  frequency  wave  will  couple  to  a  counterpropagating 
negative- frequency  wave,  due  to  the  nonlinear  interaction,  and  this  wave  exits 
the  medium  at  the  boundary  z  -  -A.  We  call  this  the  nonlinear  wave  (ni-wave), 
and  its  wave  vector  k^^  is  approximately  equal  to  k^.  Taking  into  account  the 
frequency  mismatch,  we  find  for  the  z-component  of  k^^ 
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,  ,'272  72 

k  .  -  Jp  K  -  k||  -  -k 

ni,z  II  pc,z 


(4.13) 


The  total  positive- frequency  field  in  z  <  -A  then  becomes 


E(r.t)<+)  -  E. 

-  -  ka  ko-ka 


+  E,  N*  e,  e 
kff  k^T-k  .a 
-  -  -ni 


(4.14) 


in  terms  of  the  Fresnel  transmission  coefficients  T,  and  N,  for  the  t-wave 

ko  kc7 

and  ni-wave,  respectively.  Notice  that  we  label  all  Fresnel  coefficients  with 
the  wave  vector  k  of  the  incident  wave,  rather  than  with  the  wave  vector  of 
the  corresponding  wave.  This  will  turn  out  to  be  most  practical  later  on. 

For  a  given  incident  plane  wave  we  find  a  total  of  five  plane  waves 
outside  the  medium:  three  in  z  >  0  and  two  in  z  <  -A.  It  is  easy  to  show 
chat  this  is  the  most  general  solution  of  Maxwell's  equations  (or  Heisenberg's 
equations  for  quantum  fields) ,  due  to  the  restrictions  put  on  the  wave  numbers 
by  the  dispersion  relation  in  vacuum,  by  the  causality  requirement  that  all 
waves  must  travel  away  from  the  medium,  and  by  the  phase -matching  conditions 
in  z  -  0  and  z  -  -A.  In  the  same  way  it  can  be  shown  that  the  most  general 
solution  for  the  waves  inside  the  medium  consists  of  four  sets  of  two 
counterpropagating  waves.  This  brings  the  total  count  of  the  number  of  waves 
to  thirteen,  which  are  all  present  due  to  the  excitation  by  a  single  incident 
plane  wave.  By  matching  these  waves  across  the  boundaries  at  z  -  0  and  z  - 
-A,  we  can  derive  explicit  expressions  for  the  Fresnel  coefficients  [24].  In 
this  paper  we  shall  consider  these  coefficients  as  given  parameters  of  the 
problem.  Figure  1  illustrates  the  situation. 
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V.  VACUUM  FIELD 

With  the  plane-wave  solution  from  the  previous  section,  we  can  now 
construct  the  explicit  form  of  the  vacuum  field  E^(r,t),  and  because  the 
dipole  which  couples  to  this  field  is  situated  in  2  >  0,  we  shall  consider 
in  this  region  only.  To  this  end,  we  first  recall  that  the  positive -frequency 
part  of  the  electric  -  field  operator  in  empty  space  is  given  by 


^  /-i  _  _ 


-OJt) 


VxJ 


(5.1) 


with  (ji  “  ck,  V  the  quantization  volume,  and  a^  the  annihilation  operator  for 

a  photon  in  the  mode  ko.  For  the  unit  polarization  vectors  e.  we  take  e. 

-  -ka  -ks 

from  Eqs .  (4.4)  and  (4.5),  respectively.  Then  we  introduce  the 


and  e, 
-kp 

operators 


o  - 


(5.2) 


which  gives 


ko 


-wt) 


(5.3) 


This  shows  that  E  (r,t)^^^  is  a  summation  over  ka  of  plane  waves  which  have 
— e  —  — 

exactly  the  form  (4.1)  of  an  incident  monochromatic  plane  wave.  Since  the 
plane-wave  solution  from  the  previous  section  is  the  only  one  admitted  by  the 
Maxwell -Heisenberg  equations,  we  conclude  that  every  ka-term  in  Eq.  (5.3)  must 
give  rise  to  an  excitation  of  a  fundamental  plane -wave  mode. 
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In  Sec.  IV  we  considered  a  plane  wave  which  was  incident  from  the  region 
z  >  0,  corresponding  to  a  negative  value  of  k^.  In  the  summation  (5.3),  half 
the  number  of  wave  vectors  k  has  a  negative  z -component,  and  for  these  modes 
we  can  simply  take  Eq .  (4.10),  summed  over  Vjj  ,  for  the  field  In  z  >  0.  The 
other  half  of  the  terms  in  Eq .  (5.3)  have  k^  >  0 ,  and  they  correspond  to 
incident  waves  from  the  region  z  <  -A.  In  z  >  0,  these  waves  generate  t-waves 
and  ni -waves  in  analogy  to  Eq.  (4.14).  Therefore,  the  vacuum  field  in  z  >  0 
is 


.  ,  i(k  ‘r-wt) 

Mr.t)  -  )  ^  ,  e  ^  '  , 

V.a  ~ 


-Y 


\  ^  ~k  a 


-i(k  T+pwt) 
~pc  - 


■^2  ® 
Vo 


i(k*r-wt) 


F  2k  ,, 

f-  -  -  -ni 

Vxj 


+  H.c. 


(5.4) 


where  a  prime  (double  prime)  on  a  summation  sign  indicates  that  the  summation 

runs  over  k  values  with  k  <  0  (k  >0)  only.  The  summations  on  the  right- 
-  z  z 

hand  side  of  Eq.  (5.4)  form  together  the  positive- frequency  part  of  the  field, 

and  the  addition  of  the  Hermitian  conjugate  then  yields  the  total  vacuum 

field.  The  primes  on  T,'  and  N,’  are  a  reminder  that  these  are  the  Fresnel 

ka  lu7 
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coefficients  for  a  wave  which  is  incident  from  the  region  z  <  -A,  and  is 
identical  in  form  to  expression  (4.1).  The  absolute  value  of  T  and  T'  is  of 
course  the  same  when  both  incident  waves  have  the  same  freuquency,  angle  of 
incidence,  and  polarization,  but  their  phase  will  be  different.  For  a 
transparent  medium  (unit  dielectric  constant)  the  r-waves  and  ni - 
waves  disappear  to  a  good  approximation,  which  gives  ® 

(5.4).  The  t-waves  in  z  >  0  remain  present  in  this  limit,  and  they  were 
called  'quantum  noise'  by  Gaeta  and  Boyd  [28],  following  Caves  [29].  For  a 
transparent  medium,  and  without  the  nonlinear  interaction,  we  obviously  have 
1  (because  then  E  —  E  ) .  It  can  be  shown  [24,30]  that  the  nonlinear 
interaction  always  gives  ^  provided  that  the  dielectric  constant  is 

unity . 

VI.  VACUUM  CORRELATION  FUNCTIONS 

With  the  explicit  expression  (5.4)  for  the  vacuum  field,  we  can  now 
evalute  the  vacuum  correlation  functions  from  Eq.  (3.6).  We  set  r  -  h, 

and  t  ••  T  and  t  —  0  for  the  i  and  the  j  component  of  E  ,  respectively. 
Multiplying  both  expressions  and  taking  the  vacuum  expectation  value  then 
yields 
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r-,  i(k  •h+pwr) 

-)  JTW'k.  Sk  ,  '  ■ 

,  o  -  -pc 
ka 


+ 


ik«h 
e  -  -  + 


ik  ‘h 

-r  -iwr 
e  )e 


’i 


-ik  ‘h  ,  -ik  ‘h 

X  (P*  e  e  +  (e,  e  +  RJ  e  e  ^  )). 

lo7  -k  £7  -k£7  Tc<7  -k  a  j 

-pc  -  -  -r 


ka 


i(k  -h+pwr)  . /I  u  \ 

cj  ...  ~ni  -  L(k*h-<JT), 

T — —  { N,'  e  e  -I-  e,  T,'  e -  )  . 

2e  W  kc7  -k  .£7  -k£7  k£7  i 

o  -  ~ni  -  - 


{N/*  e. 

k£7  -k  .£7 

-  -ni 


-ik  ‘h 
-ni  - 


~k£7 


T'  * 
k£7 


•  ik*h j 


(6.1) 


Next  we  have  to  substitute  the  values  of  the  various  wave  vectors,  as  they 

were  introduced  in  Sec.  IV.  We  notice  that  all  wave  vectors  in  the  exponents 

appear  in  the  combination  k.«h  -  hk.  ,  so  that  we  only  need  their  z- 

-1-1, z 

components.  At  this  stage  we  make  a  slight  approximation.  The  Fresnel 

coefficients  P,  and  N,'  only  deviate  considerably  from  zero  when  the 

frequency  w  of  the  incident  wave  is  in  close  resonance  with  the  pump  frequency 

w.  Therefore,  in  factors  multiplying  these  Fresnel  coefficients  we  can  safely 

replace  w  by  w .  This  amounts  to  the  approximations  k  =•  k  ,  k  .  “ 

f  Zt  Z  XiA  f  z 

-k  ,  e,  “  e,  and  e,  =•  e,  ,  but  the  full  resonant  behavior  of  the 
z  -k  a  -ka  -k  .a  -k  a 

-pc  -  -ni  -r 

Fresnel  coefficients  as  a  function  of  cj  will  be  retained.  Then  we  obtain 
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-  2tk*h,  -  iu)f 


c  /  \  \  ^  XpCi)T  ,  .T>  xwr  . 

“  Z  2e  XV  Ska  "  (Ska  ^  \a  Sk_a  ®  ’i 


■ij 


ka 


--r 


-r 


'-^  -2ik*h  ipwT  _,  -iwr, 

^  I  iT^  (^ka  Sk  a  ^  ~  ^  ^  ^ka  ^ka  ®  'i 

f-  o'  -  -r  -  - 

ka 


, .  L  2 ik"h  „ .  . 

X  N,'  *  e,  e - +  T,'  *  e  )  . 

ka  -k  a  ka  -ka  i 

-  -r  -  -  -' 


(6.2) 


Notice  that  we  have  not  set  p  -  1  in  exp(ipwr),  since  this  would  affect  the  r- 
dependence  of  f_(r). 

As  the  next  step  we  want  to  perform  the  summations  over  k,  with  the  help 
of  the  identity 


I  d\  (...) 

k  <0 


(6.3) 


and  similarly  for  the  double-prime  summation.  The  integral  is  most  easily 
carried  out  in  spherical  coordinates,  for  which  we  have 


(6.4) 


and  for  the  summation  over  k^  >  0  we  replace  the  integration  limits  (x/2,x) 
for  the  integral  by  (0,x/2)  For  the  k-vector  we  have 
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k  -  ke 
-r 


(6.5) 


and  with  Eqs.  (4.4)  and  (4.5)  we  find  for  the  polarization  vectors 


e,  -  e  -  -e 
-ks  -k  s  ~<p 
~  ~r 


(6.6) 


e,  “0., 

-kp  ~d 


(6.7) 


e,  -  -e,  -2sin^e 
-k  p  ~d  ~z 

-r 


(6.8) 


in  terms  of  the  standard  spherical  unit  vectors  e  ,  e  and  e  . 

— r  “d 

When  we  make  the  substitutions  (6.4)-(6.8)  into  Eq.  (6.2)  then  we  get  a 
very  lengthy  expression  for  fj^j(r),  and  for  every  one  of  the  nine  combinations 
(i,j)  this  expression  is  different.  A  great  simplification  arises  from  the 
fact  that  the  Fresnel  coefficients  only  depend  on  k  and  6  (or  the  frequency  of 
the  incident  wave  and  the  angle  of  incidence),  but  not  on  Therefore,  the 
only  ^-dependence  of  the  integrand  enters  through  the  polarization  vectors 
from  Eqs.  (6.6)-(6.8),  and  the  integration  over  (f>  can  be  performed  immediately 
in  each  of  the  nine  combinations.  We  then  find  that  this  integral  equals  zero 
whenever  i  j  ,  and  that  i  -  j  -  x  gives  the  same  result  as  i  -  j  -  y. 
Consequently,  we  can  write 


-  ^ij  ^ii^"^  ’ 


(6.9) 


and  for  the  three  nonzero  correlation  functions  we  have 
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fl|(»-)  -  f  (r)  -  f  (r)  ,  (6.10) 

II  XX  yy 

fl(0  -  f^^(r)  ,  (6.11) 

which  defines  fx(r).  It  can  be  shown  [31]  that  results  (6.9)-(6.11) 

are  not  just  a  lucky  simplification  for  the  present  problem,  but  that  these 
relations  can  be  derived  in  general  from  the  invariance  of  the  vacuum  for 
rotations  about  the  z-axis,  and  for  reflections  in  a  plane  through  the  z-axls. 
For  a  spherically- symmetric  vacuum  (like  empty  space),  we  find  in  addition 
that  fj^(r)  must  necessarily  equal  f||(r). 

For  the  integration  over  6  we  change  the  integration  variable  to  the 
angle  of  incidence  b  for  a  given  value  of  k,  and  the  relation  between  b  ^  and 
b  is 

IT  -  5  ,  for  k  <0 

z 

5.  -  ■  (6.12) 

b  .  for  k  >  0 

Then  we  know  that  Fresnel  coefficients  only  depend  on  b^  through  costf^  (or 
sin^^) ,  because  their  angle-of-incidence  dependence  is  purely  geometrical  and 
can  be  expressed  in  cross  and  dot  products  between  the  various  wave  vectors, 
polarization  vectors  and  e^ .  We  write 

P,  -  P  (cj,  costf  .  ) 

Va  a  I 


(6.13) 


22 


with  w  >  0  the  frequency  of  the  incident  wave,  and  similarly  for  the  other 
Fresnel  coefficients.  For  the  k-  and  ^ - i ntegration  we  then  change 
integration  variables  according  to  w  -  ck  and  u  -  cosfl^,  respectively.  For 
the  two  vacuum  correlation  functions  we  then  obtain 
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where  all  Fresnel  coefficients  have  (w,u)  as  their  arguments.  These 
coefficients  depend  in  complicated  way  on  w  and  u,  in  general,  and  it  does  not 
appear  possible  to  evaluate  the  integrals  in  Eqs.  (6.14)  and  (6.15). 


VII.  LAPLACE  TRANSFORM 

The  vacuum  correlation  functions  f.j(r)  only  enter  the  spontaneous -decay 
operator  through  from  Eq.  (3.5).  When  we  expand  exp(iL^r)  onto  the  energy 
eigenstates  of  the  atom,  then  this  factor  effectively  turns  into  exp(iwr)  with 
u)  the  energy  separation  between  two  levels,  divided  by  h-  Consequently,  the 
vacuum  correlation  functions  only  appear  as  their  Laplace  transform 


f  (w) 


iur 

e 


(7.1) 


where  u  can  be  either  positive  or  negative.  Since  the  factors  in  Eqs.  (6.14) 
and  (6.15)  which  depend  on  r  are  all  exponentials,  we  can  evaluate  the  Laplace 
transforms  with  the  identity 


f*  dr  -  irS(<^±u’)  +  P 


(7.2) 


where  P  stands  for  principle  value.  We  shall  omit  the  principle -value  part, 
since  it  leads  to  a  small  level  shift  and  not  to  spontaneous  decay.  After 
integration  over  r  we  can  also  perform  the  w- integration,  which  finally  gives 

r;(w)(x  (w)  +  z  (w) )  ,  cj  >  0 

a  a 

(7.3) 

f;(2w-Ki)){y  (2uH-u>)  +  z  (2w-(-w)*)  ,  w  <  0 

a  Q 


f  (w)  -  i 
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with  a  -  Here  we  introduced 


rj(w)  - 


(7.4) 


Sire  ^c' 
o 


and  the  six  parameter  functions 


^±(<^)  “  4  +  |T'l 


Mu|2  ^  ,.,,,2, 


(7.5) 


Xii  (ta) 


[  du  (  I  1  +  R  +  u^i 

J  n  ^ 


-  k  e 
P 


i^U|  2 


+  1t;|^  ^  u^lT’l^l 


yi(a,)  -  J  f  du  (l-u2)(|P  +  iN'I^l 

Jq  p  p 


yii(<^) 


-  t  du  {  |pj2  +  |n;I^  +  VISIN' 1^)  . 


ZiCw)  -  7  f  du  (l-u^)(P*  +  e^^^(P*R  +  T'N'*))  , 

^Jn  P  PPPP 


(7.6) 


(7.7) 


(7.8) 


(7.9) 


Z||(u;)  -  I  I  du  {P*  +  u^P*  +  e’’^“(P*R 
II  3  J-  s  p  ^  s  s 


u  P*R  +  T'N'* 
P  P  s  s 


u  T'N'*))  . 

P  P 


(7.10) 


which  are  defined  for  w  >  0  only.  All  Fresnel  coefficients  have  (ca.u)  as 
their  argument,  and  the  parameter  ^  is  defined  by 
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0 


2a)h 

c 


(7.11) 


Notice  that  the  atom-surface  distance  h  only  enters  via  the  parameter  and 

therefore  the  two  parameter  functions  yx(w)  and  y||  (w)  are  independent  of  this 

distance.  Furthermore,  the  functions  x  (w)  and  y  (w)  are  real  and  non- 

a  ■'a 

negative,  but  z^((j)  and  Z||((j)  are  complex-valued,  in  general. 


VIII.  PARAMETER  FUNCTIONS 

Spontaneous  decay  of  any  atom  near  any  linear  or  nonlinear  medium 
occupying  the  region  0  >  z  >  -A  is  determined  completely  by  the  two  functions 

fj^(w)  and  fi|(w),  which  can  be  parametrized  by  the  six  functions  x  (w)  ,  y  (w)  , 

''  % 

z^(w).  These  six  parameter  functions  are  again  determined  by  the  eight 
Fresnel  coefficients  R  ,  P  ,  T' ,  N'  (c?  -  s,p),  and  the  atom-surface  distance 
parameter  As  soon  as  the  Fresnel  coefficients  for  a  certain  configuration 
are  given  as  a  function  of  the  cosine  of  the  angle  of  incidence  u,  then  the 
integrals  in  Eqs.  (7. 5) -(7. 10)  can  be  evaluated,  in  principle.  The  dependence 
on  the  frequency  w  is  only  parametric,  where  £*>  equals  a  Bohr  frequency  of  the 
atom.  In  order  to  illuminate  the  significance  of  the  parameter  functions,  we 
work  out  some  important  limiting  cases. 


A.  Empty  space 

Without  a  medium  we  simply  have  an  atom  in  r  -  h ,  and  surrounded  by 
empty  space.  Then  the  Fresnel  coefficients  obviously  are 

R-P-N'-O  ,  T'-l  ,  (8.1) 

a  o  o  o 


and  the  parameter  functions  are  readily  found  to  be 
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x-1  ,  y-z-0 

a  a  a 


(8.2) 


They  are  the  same  for  q  -  and  a  -  || ,  which  indicates  the  spherical  symmetry 
of  the  vacuum,  and  they  are  independent  of  u  (and  of  h) . 


B.  Linear  medium 

Wnen  the  medium  is  a  dielectric,  a  metal,  a  thin  film  or  a  substrate,  or 
any  composition  for  which  reflection  and  refraction  is  induced  by  the  first- 
order  susceptibility  only,  then  every  plane  wave  has  the  same  frequency  as  an 
incident  plane  wave,  and  consequently  we  have 

P-N'-O,  Rj'O,  T'^O,  (8.3) 

a  a  a  a 

in  general.  Then  we  find  immediately 


X  0 
a 


y  -  z  -  0 
■'a  a 


and  the  vacuum  correlation  functions  become 


(8.4) 


fjcj) 

a 


rf(u>)x  (w) 
a 


0 


w  >  0 


w  <  0 


(8.5) 


which  will  be  different  for  the  two  values  of  a.  For  positive  frequencies, 
f^(w)  is  real  and  positive,  but  for  negative  frequencies  f^(w)  vanishes 
identically.  Conversely,  f^(<^)  can  only  acquire  a  finite  value  for  w  <  0  due 
to  a  nonlinear  process  in  the  medium. 
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C.  Perfect  conductor 

The  most  simple  and  nontrivial  linear  medium  is  the  perfect  conductor 
(mirror) ,  for  which  we  have 

T'-O,  R--1,  R-1,  (8.6) 

asp 

in  addition  to  Eq .  (8.3).  Then  the  integrals  over  u  can  be  found 

analytically,  and  we  get 

-  1  -  3(^  -  ,  (8.7) 

X|,<»,  -  1  .  ^  .  (8.8) 

which  is  a  well-known  result.  It  exhibits  the  famous  oscillatory  structure  as 
a  function  of  h  (or  0)  ,  and  in  the  limit  h  -♦  ®  we  find  Xj^  -  X||  -  1.  This  is 
the  empty- space  limit  from  subsection  A,  and  it  displays  that  far  away  from 
the  surface  the  spontaneous  decay  is  not  affected  anymore  by  the  presence  of 
the  medium. 


D.  Transparent  PC 

For  most  four-wave  mixing  PC's  one  uses  a  transparent  crystal  (unit 
dielectric  constant) ,  which  has  the  advantage  that  the  specular  wave  is 
(almost)  absent  and  that  the  phase-conjugate  reflectivity  can  be  high.  In 
principle,  there  is  no  limit  on  the  value  of  P  ,  i.e., 
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0  <  ®  .  (8.9) 

and  for  the  specular  waves  we  have 

R  -  N'  -  0  .  (8.10) 

a  a 

The  t-waves  are  always  there,  and  their  Fresnel  coefficients  are  related  to  P 

o 

according  to 

=  1  .  (8.11) 

where  both  T  and  P  pertain  to  the  same  incident  plane  wave.  It  can  be 
shown  [24]  that  relation  (8.11)  holds  for  any  angle  of  incidence  and  any 
frequency  mismatch  between  w  and  w.  In  particular,  Eq.  (8.11)  implies 

>  1  .  (8.12) 

indicating  that  the  transmitted  wave  is  always  amplified  in  the  four-wave 

mixer,  as  compared  to  the  incident  wave.  Equation  (8.9)  allows  for  values  of 
2 

I  I  larger  than  unity,  in  which  case  also  the  pc-wave  is  stronger  than  the 
incident  wave.  Amplification  in  a  four-wave  mixing  PC  has  been  observed 
experimentally  indeed  [32,33].  It  is  also  worth  noting  that  relation(8 . 11)  is 
responsible  for  the  generation  of  two-photon  coherent  states  in  four-wave 
mixing  [34,35]  . 

For  this  most  important  and  realistic  example  we  find  that  x^  and  y^  are 
related  as 


I 


I 
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X  (w)  -  1  +  y  (oj) 
a  a 


(8.13) 


with 


y_L(a,)  -  J  [  du  (1  -  u2)|P  |2 


(8.14) 


^11 


(w) 


-f  du  I|pj2  *  „2|p^|2|  , 


(8.15) 


2  2 

where  we  have  used  Eq .  (8.11)  to  eliminate  |t  I  in  favor  of  |P  I  .  It  should 

'a'  '  £7 ' 


be  noted  that  T'  from  Eqs .  (7.5)  and  (7.6)  is  not  equal  to  T  in  Eq.  (8.11), 


although  their  absolute  values  are  equal  when  they  both  pertain  to  the  same 


angle  of  incidence.  For  z  we  find 

a 


zi(w) 


-If' 

"Jo 


du  (1  -  u  )P*  , 

P 


(8.16) 


Z||  (w) 


-  c 

«  Jo 


du  { P*  +  u  P* ) 
s  p 


(8.17) 


In  the  expression  for  f^(w),  Eq.  (7.3),  we  always  get  y^  and  in  the 

combination  y^  +  z^.  From  Eqs.  (8.14)  and  (8.15)  we  see  that  >  0,  and  that 

also  the  integrand  is  non-negative  for  every  value  of  u.  Therefore,  the 

contributions  from  all  plane  waves  with  different  angles  of  incidence  add  up 

constructively.  This  in  contrast  to  Eqs.  (8.16)  and  (8.17)  where  the  values 

of  P*  appear,  rather  than  their  absolute  values.  The  Fresnel  coefficients  P 
a  a 

vary  rapidly  as  a  function  of  the  angle  of  incidence  [24]  ,  and  in  an 
oscillatory  fashion  around  P^  -  0.  Integrated  over  the  angle  of  incidence,  as 
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in  Eqs.  (8.16)  and  (8.17),  this  will  yield  a  very  small  value  as  compared  to 

2 

Che  integrated  intensities  |P^|  .  A  safe  approximation  is  therefore 


z  («)  “  0  . 

a 

when  used  in  Eo.  (7.3) 


(8.18) 


E.  Ideal  PC 

2  2 

When  the  values  of  |  |  and  |P  |  would  be  independent  of  the  angle  of 

2 

incidence  (and  therefore  equal  to  |P  (tj,u  -1)1  ,  because  |P  I  -  |P  I  for  u  - 

’  s  '  s  p ' 

1),  then  we  could  take  them  outside  the  integrals  in  Eqs.  (8.14)  and  (8.15), 
which  would  give 


yi(w)  -  y||(w)  -  2 


(8.19) 


For  illustrative  purposes  later  on  we  define  an  ideal  PC  as  a  transparent  PC 
for  which,  in  addition,  Eq.  (8.19)  holds.  The  reservoir  .'orrelation  function 
then  equals 


rj(u)  { 1  +  h 


w  >  0 


f  (w) 

a 


hf7(  -cj)  I  P|  ^ 


w  <  0 


(8.20) 


both  for  a  -  and  a  -  ||  .  Furthermore  we  have  assumed  that  |P^|  is  symmetric 
around  u>,  as  a  function  of  w,  which  is  quite  accurate  in  general.  We  like  to 
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emphasize  that  Eq.  (8.19)  is  not  necessarily  a  good  approximation,  and  also 

2 

that  |p|  on  the  right-hand  side  should  not  be  interpreted  as  the  reflectivity 

2  2 

at  normal  incidence.  The  value  |p1  can  be  viewed  as  an  average  of  |P^|  over 

the  angle  of  incidence,  where  the  'average'  is  defined  by  Eqs.  (8.14)  and 

(8.15).  The  approximation  is  then  that  Eqs.  (8.14)  and  (8.15)  give  comparable 

results.  In  an  approximation  where  P^(w,u)  is  assumed  to  be  independent  of  u, 

we  would  find  in  addition  that  z^(co)  -  hP* ,  according  to  Eqs.  (8.16)  and 

(8.17).  It  is  easy  to  show  that  this  leads  to  inconsistencies,  and  that  the 

appropriate  approximation  is  given  by  Eq .  (8.18).  Also  note  that  f^(<*>)  is 

identical  in  form  to  the  correlation  function  of  a  thermal  photon  reservoir 

2 

with  an  average  photon  number  of  (p|  /2. 

IX.  TWO  LEVELS 

With  the  explicit  result  for  the  vacuum  correlation  functions  from  Sec. 

VII,  we  can  evaluate  the  spontaneous -decay  operator  F  from  Eq.  (3.4)  for  any 

configuration  of  atomic  levels.  In  this  section  we  work  out  the  general  form 

of  r  for  the  case  of  two,  possibly  degenerate  or  nearly-degenerate ,  levels. 

With  hw  and  hw  as  the  energies  of  the  excited  state  and  ground  state, 

®  8 

respectively,  the  atomic  Hamiltonian  reads 

H  -  Hw  P  +  hw  P  ,  (9.1) 

a  '  e  e  g  g 

in  terms  of  the  proiectors  P  and  P  on  both  states.  The  corresponding 

eg 

Llouvillian  is 

L  a  -  [P  ,a\ 
a  o  e 


(9.2) 
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with  w  -  w  -  w  as  the  transition  frequency,  and  where  we  used  the  closure 
o  e  g 

relation 


P  +  P  -  1  (9.3) 

eg 

For  the  calculation  of  from  Eq .  (3.5)  we  need  the  expansion  of  exp(-iL  r) 
in  projectors 


-iL  r 
a 

e  a 


(j)  )r 
a 


P  aP 
a 


0 


(9.4) 


which  can  be  proven  easily  with  -  P^,  o  -  e,g.  This  exponential  acts  on 
the  Cartesian  component  of  the  dipole  operator,  which  gives 


-  iL  r  -ia>  r 

!  ^  p.  -  e  °  ?  ^i.P 

i  e^i  g 


i(<>  T 

e  °  P  /i.P 
g  1  e 


(9.5) 


and  where  we  used 


(9.6) 


This  relation  follows  from  the  fact  that  an  atom  cannot  have  a  permanent 
dipole  moment.  Substitution  of  the  right-hand  side  of  Eq .  (9.5)  into  Eq . 
(3.5)  then  yields 


Q,  -  f,,(-w  )P  M.P  +  )P  M,P 

i  ii  o  eig  ii  o  gie 


(9.7) 


For  molecular  transitions,  Eq.  (9.6)  does  not  hold  since  a  molecule  can  have  a 


permanent  dipole  moment.  Nevertheless,  the  expression  for  remains  the  same 
because  the  terms  with  become  proportional  to  f^j^(O),  which  equals 

zero . 

Then  we  substitute  and  its  Hermitian  conjugate  into  Eq.  (3.4),  work 
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out  the  commutator,  and  drop  the  nonsecular  terms,  which  gives  for  F 


Fcr  -  )  lM?[P  f..(-w  )  +  P  f..(w  )]a 
/_  L  g  11  o  e  11  o 


+  a[?  f..(-(.j  )*  +  P  f..(aj  )*1m- 
g  11  o  e  11  o  1 


-  2/1.  [P  a?  Re  f..(-a)  )  +  P  oP  Re  f .  .  (uj  )]/i.) 
1  g  g  11  o  e  e  11  o  -^1 


(9.8) 


Here  we  have  used  identities  like 


P  /i.P  - 
e'^l  g 


P  /i  .  -  /i  .  P 

e^l  ^1  g 


(9.9) 


which  easily  follow  from  Eqs.  (9.3)  and  (9.6). 

In  Eq.  (9.8)  we  have  expressed  F  In  terms  of  the  Cartesian  components 
of  the  dipole  operator  and  the  projectors  P^  and  P^  on  both  levels,  which  is 
the  most  compact  representation  of  F.  For  many  purposes,  however,  it  is  more 
convenient  to  express  F  in  atomic  lowering  and  raising  operators,  rather  than 
in  projectors.  The  lowering  and  raising  parts  of  /i  are 
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-  P  mP 

-  6=  g 


respectively,  and  the  total  dipole  moment  is 


(9.11) 


(9.12) 


Then  we  substitute  its  i-th  component  into  Eq .  (9.8)  and  use  the  notation  with 
-1  and  11  ,  which  gives 


To  -  I  f  i|  ( -oi  )/i  II 

1  o  - 


(  +  ) 


+• 


f  II 

11  o  - 


■t-  fi( 'PX  ^  +■  P±(‘^o)Mi 


+  o(  f  II  ( -w^)*P|j  ^*4|  ^  ^*-|j 

+  fl( *^1  ^  +  Pl(“Q)*ifi 

-  2  ^  fj|(-u.^)  -  2  ^  pJ^^p[‘^Re  fj|(u>^) 

i-x ,  y  i-x ,  y 

-  Re  f,(-w  )  -  f.(w  )  .  (9.13) 

z  o  z  z  -*-  o 


We  notice  that  every  term  is  proportional  to  from  Eq .  (7.4),  but  since 
there  is  only  one  transition  frequency  for  a  two-level  system  we  simply  write 
V  ■ 


35 


X.  TWO -STATE  ATOM  AND  IDEAL  PC 

In  order  to  discover  the  physical  mechanisms  of  spontaneous  decay  near  a 
PC,  and  their  relation  to  fluorescent  emission,  we  consider  the  model  case  of 
a  two-state  atom,  in  combination  with  the  ideal  PC  for  which  the  vacuum 
correlation  function  is  given  by  Eq .  (8.20).  With  |  e>  and  |  g>  as  the  wave 
functions  of  the  excited  state  and  ground  state,  respectively,  the  projectors 
on  these  states  are 


Then  expression  (9.13)  for  F  can  be  simplified  considerably,  and  in  its  most 
transparent  form  it  reads 
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fa  -  4  A  {P  t7  +  aP  -  2d ‘ad) 
e  e  e 


+  4  A  { P  a  +  ctP  -  2cL7d  ) 
g  g  g 


(10.5) 


Here  we  have  introduced  the  real  and  positive  parameters 


A^  -  A(1  +  hlPp)  , 


(10.6) 


A  -  h  aIpI' 
g 


(10.7) 


whose  significance  will  become  clear  in  due  course. 

The  equation  of  motion  for  the  atomic  density  operator  is  given  by  Eq. 
(2.6)  (with  p  -  p)  .  Most  interesting  is  the  time  evolution  of  the 
populations  of  the  two  levels,  defined  by 


n  (t)  -  Tr  P  p(t) 
a  a 


a  -  e,g  . 


(10.8) 


When  we  multiply  Eq .  (2.6)  on  the  left  by  P^,  take  the  trace,  use 


TrP(Lcr)-0  ,  a-e,g 

Cl  3 


(10.9) 


and  work  out  the  result  for  TrP^(ra) ,  as  it  follows  from  Eq.  (10.5),  then  we 


obtain 


--An  +  A  n 
e  e  g  g 


(10.10) 
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dn 
— g 
dt 


-An  +  A  n 
g  g  e  e 


(10.11) 


First  we  notice  that  the  right-hand  sides  of  these  equations  differ  only  by  a 
minus  sign,  and  therefore  n^  +  is  independent  of  time.  The  normalization 
is 


n  (t)  +  n  (t)  -  1  .  (10.12) 

e  g 

which  is  nothing  but  Trp(t)  -  1.  Equations  (10.10)  and  (10.11)  relate 

populations  only,  and  they  have  the  form  of  a  set  of  loss/gain  equations.  The 

physical  interpretation  of  these  rate  equations  is  that  equals  the  number 

of  le>  -»  I  g>  transitions  per  unit  of  time,  and  An  is  the  number  of  |  g>  le> 

g  g 

transitions  per  unit  of  time.  Therefore,  A  and  A  are  the  rate  constants  for 

e  g 

|e>  -►  I  g>  and  |  g>  -•  |  e>  transitions,  respectively,  and  they  describe 

spontaneous  transitions  by  this  atom  in  the  vacuum  near  a  PC.  In  the  limit 
2 

|p|  -►  0  the  value  of  A  reduces  to  A,  whereas  A  vanishes.  Then  we  have 

e  g 

ordinary  spontaneous  decay,  with  A  as  the  rate  constant.  Due  to  the  presence 

of  the  PC,  the  value  of  A  increases  to  A  ,  which  shows  that  the  PC  has  the 

e 

effect  of  enhancing  the  spontaneous -decay  rate.  More  interesting  is  that  A 

g 

acquires  a  finite  value,  which  implies  that  there  are  spontaneous  excitations 

of  the  atom  in  vacuum,  at  a  rate  An.  It  is  also  worth  noting  that  the 

g  g 

Increase  of  both  A  and  A  ,  as  compared  to  their  empty- space  values,  equals 
6  g 

2 

hA  I  P  I  This  situation  is  reminiscent  of  the  coupling  of  a  small  system  to  a 
nonzero- temperature  reservoir,  where  the  increase  of  the  relaxation  constants 
due  to  the  finite  temperature  is  also  equal  for  decay  and  excitation.  An 
obvious  interpretation  of  Eqs .  (10.6)  and  (10.7)  is  that  the  atom  will  decay 
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spontaneously  with  a  rate  constant  equal  to  A,  but  that  in  addition  we  have 
stimulated  |e>  -►  |  g>  and  |  g>  |  e>  transitions,  which  are  induced  by  the  PC, 
and  which  increase  the  rate  constants  for  both  down  and  up  transitions  by  the 
same  amount. 


z 


XI.  COMPARISON  WITH  FLUORESCENT  EMISSION 

As  mentioned  in  the  Introduction,  spontaneous  decay  is  intimately 
related  to  spontaneous  emission  because  of  energy  conservation.  For  the 
present  situation,  however,  this  relation  is  not  too  obvious  as  follows  from 
the  fact  that  there  is  the  possibility  of  spontaneous  excitation.  A  | g>  -»  |e> 
transition  requires  the  supply  of  an  amount  of  energy  equal  to  .  and  the 
process  is  of  course  the  absorption  of  a  photon.  Only  the  PC,  with  its  two 
pump  beams  can  provide  this  energy. 

Elsewhere  [19]  we  have  derived  an  expression  for  the  emitted  fluorescent 
power  dW/dt  by  an  independent  method  (solving  Maxwell's  equations  without  an 
explicit  quantization  of  the  radiation  field).  We  only  considered  the 
emission  into  the  positive  z-direction,  since  only  this  field  can  be  measured 
directly  by  a  photomultiplier  in  the  far  zone.  For  the  number  of  emitted 
photons  per  unit  of  time  into  the  z  >  0  direction  we  found 


1  ^ 

hu)  dt 
o 


HAn  +  A  n 
e  g  g 


(11.1) 


2 

with  A^  -  hA|p|  That  the  situation  is  more  complicated  than  for  ordinary 
fluorescence  is  illustrated  with  the  four  diagrams  of  Fig.  2.  The  spontaneous 
contribution  A  to  A^ ,  Eq.  (10.6),  corresponds  to  ordinary  spontaneous  decay, 
accompanied  by  the  emission  of  a  fluorescent  photon.  This  process  happens  at 


I 
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a  rate  An^,  and  every  transition  |  e>  -♦  |  g>  produces  a  photon.  However,  as 

shown  in  diagrams  (1)  and  (2),  only  half  the  number  of  photons  travels  towards 

the  detector  in  z  >  0 ,  and  this  gives  the  term  ^4An^  on  the  right-hand  side  of 

Eq .  (11.1).  Then,  the  stimulated  contribution  to  the  relaxation  constants 

equals  hA | P |  for  both  up  and  down  transitions.  Diagram  (3)  of  Fig.  2 

illustrates  the  stimulated  | e>  -►  |g>  transition  and  the  corresponding  photon 

emission.  The  rate  of  this  process  is  proportional  to  n^,  but  the  n^-term  in 

Eq.  (11.1)  is  already  accounted  for.  Since  energy  conservation  requires  the 

emission  of  a  photon,  we  conclude  that  this  photon  must  have  frequency  «,  and 

is  therefore  indistinguishable  from  the  photons  already  present  in  the 

vacuum.  This  situation  is  analogous  to  atomic  decay  in  a  f  inite -  temperature 

vacuum,  where  the  emitted  photons  in  a  stimulated  transition  are 

indistinguishable  from  the  black-body  radiation.  The  stimulated  |  g>  -♦  |  e> 

excitations  are  shown  in  diagram  (4),  and  they  correspond  to  the  term  An  - 

&  S 

2 

hA|p|  n  on  the  right-hand  side  of  Eq.  (11.1).  Since  this  term  is 
& 

proportional  to  n  ,  the  atom  must  be  originally  in  its  ground  state.  An 

6 

absorption  of  a  vacuum  photon  with  frequency  «  brings  the  atom  to  its  excited 
state,  and  a  subsequent  decay  produces  a  fluorescence  photon.  It  follows  from 
the  rate  equations  (10.10)  and  (10.11)  that  the  atom  must  be  in  state  | e> 
after  completion  of  the  process,  and  therefore  a  second  photon  with  frequency 
w  will  be  absorbed.  This  three-photon  process  gives  a  fluorescent  photon  with 
frequency  2a)-(a^,  whereas  the  spontaneous -decay  process  from  diagram  (1) 
produces  a  photon  with  frequency  With  frequency  selective  detection  it 

should  be  possible  to  discriminate  between  the  two  processes  in  an 
experimental  observation  of  the  fluorescence. 
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XII.  RELAXATION 

The  set  of  equations  ( 10 . 10)  ,  ( 10 . 11)  is  easily  solved.  For  a  given 
population  of  the  excited  state  at  t  -  0,  the  solution  for  t  >  0  is 


n  (t)  -  ^ 

e  A  +  A 


S 


-(A  +A  )t  A 

e  g 


(12.1) 


and  the  population  of  the  ground  state  follows  from  Eq.  (10.12).  For  t  »  0 
we  obtain  the  steady-state  value 


ti  (“)  -  I 
e  A  +  A 


g  1  +  IpI 


(12.2) 


Without  a  PC  the  long-time  limit  is  always  n^(“)  -  0,  but  due  to  stimulated 

excitations  from  the  ground  state  there  is  a  finite  probability  to  find  the 

atom  in  its  excited  state.  Notice  that  this  probability  is  determined  by  the 

phase  -  conj  ugate  reflectivity  |p|  only,  and  that  it  is  independent  of  A  (and 

thereby  independent  of  the  dipole  moment  fi  ) .  It  can  be  shown  [37],  however, 

—eg 

that  this  is  an  artefact  of  a  two-state  model.  Also  notice  that  the  time 
constant  of  the  relaxation  is  the  inverse  of 


A 

e 


+  A  -  A(1  +  |P|^) 
g 


(12.3) 


which  is  larger  than  A.  Consequently,  an  atom  near  a  PC  relaxes  faster  than 
in  empty  space. 
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In  the  steady  state  the  photon  emission  rate  from  Eq.  (11.1)  becomes 


1 


o 


®) 


-  a|p| 


2  3  +  |P|^ 

1  +  |P|2 


(12.4) 


2 

which  is  nonzero  for  |p|  r*  0.  An  atom  near  a  PC  continuously  emits 
radiation,  both  due  to  ordinary  decay  (diagram  (1)  form  Fig.  2)  and  due  to 
three -photon  processes  (diagram  (4)).  From  Eq.  (12.2)  we  derive 


A  n  (®)  - 
g  g 


(2  +  pr)^^ 


An  (» 
e 


')) 


(12.5) 


2 

which  shows  chat  the  three-photon  processes  produce  a  factor  of  2  +  |p|  more 
photons  in  z  >  0  chan  ordinary  decay  from  the  excited  state. 

The  equation  of  motion  (2.6)  is  also  easily  solved  for  the  time 
evolution  of  the  coherences.  It  appears  that  the  equation  for  <e|p(t)|g>  does 
not  couple  to  the  rate  equations  for  the  populations,  and  the  general  solution 
is  found  to  be 


<e|p(t) |g> 


e 


-h(A  +A  )c 
e  g 


-  iw  t 

e  °  <e|p(0)lg> 


(12.6) 


for  a  given  initial  state  p(0).  The  relaxation  time  for  the  coherences  Is 
twice  as  long  as  the  relaxation  time  for  the  populations,  as  could  be  expected 
from  the  analogy  with  thermal  relaxation. 


XIII.  DEGENERATE  TWO -LEVEL  ATOM  AND  TRANSPARENT  PC 

Although  in  certain  experiments  one  can  prepare  an  atom  so  as  to  behave 
effectively  as  a  two-state  atom  [38],  in  general  this  is  an  oversimplification 


A2 


of  the  situation.  In  particular,  for  an  atom  near  a  PC  we  have  the  four 

processes  of  Fig.  2  between  any  set  of  two  states  with  a  nonvanishing 

transition  dipole  moment  When  the  excited  level  is  degenerate,  for 

instance,  then  there  are  stimulated  excitations  (diagram  (4))  from  a  ground 

state  to  any  of  the  excited  states.  Also,  the  assumption  of  an  ideal  PC  is 

not  very  realistic,  although  it  is  adequate  to  study  the  fundamental  processes 

from  Fig.  2.  In  this  section  we  consider  the  realistic  situation  of  a 

transparent  PC,  as  discussed  in  Sec.  VIII. D,  and  we  assume  that  the  pump 

frequency  is  reasonably  close  to  an  atomic  resonance  so  that  yj^(w)  and 

yjj(cj)  from  Eqs .  (8.14)  and  (8.15),  respectively,  are  nonzero  for  w  -  Then 

the  vacuum  correlation  function  at  w  —  w  and  w  —  -w  is 

o  o 


f^(Wo)  -  t,(a.o)(l  +  y^(a)o))  . 


(13.1) 


with  Q  -  or  11  .  Here  it  is  also  assumed  that  |P^1  is  symmetric  around  w  in 
its  frequency  dependence.  From  now  on  we  write  rj ,  yj^  and  y|j  since  the 
frequency  is  fixed. 

With  j  and  j  the  angular  momentum  quantum  numbers  of  the  excited  level 
®  8 

and  ground  level,  respectively,  we  can  write  their 

magnetic  substates.  As  the  quantization  axis  for  the  definition  of  the 

degeneracies  we  take  the  normal  to  the  surface  (z-axis).  The  excited  level 

and  ground  level  are  (2j  +1)-  and  (2j  +  l)-fold  degenerate,  and  since  the 

®  8 

dipole  moment  matrix  element  vanishes  for  i  -  i  -  0,  the  minimum  number  of 

Jg  Jg  . 

atomic  states  is  four  (for  j  -  j  -hori  -1,  i  -Oori  -0,  j  -  1). 

e  -"g  -"e  ’  -^g  -^e  ’  ■'g  ' 
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The  generalization  of  the  atomic  raising  operator  from  Eq.  (10. 3)  for  a 
degenerate  two- level  atom  is  defined  as 


d  -  )  (imlrlim)ljm  Xj  m  I 

^  L  S  S  e  e  '-'e  e  -^g  g' 

m  m 
e  g 


r  -  -1,0.1 


in  terms  of  which 


can  be  expressed  as 


(13.3) 


M 


(-) 


/2j7^  / 


d  e* 
r-r 


(13.4) 


Here,  is  a  spherical  unit  vector  with  respect  to  the  z-axis,  and 
equals  the  Hermitian  conjugate  of  the  right-hand  side. 

The  spontaneous  -  decay  operator  F  is  expressed  in  the  Cartesian 
components  of  ^  Iri  Eq.  (9.13),  and  with  Eq .  (13.4)  and  its 
Hermitian  conjugate  we  can  simplify  the  result  considerably.  In  analogy  to 
Eq.  (10.5)  for  a  two -state  atom,  we  now  obtain 


Fa  —  h  )  A^|d  d^o  +  ad  d^  -  2d^ad  ) 

/_  T  T  T  T  T  T  T 

T 

+  4  )  A®  {d^d  a  +  ad^d  -  2d  crd^ 
/rrr  rr  rr 


(13.5) 


The  similarity  with  Eq .  (10.5)  is  even  more  obvious  if  we  write  -  dd^ , 

P  -  d^d  in  Eq .  (10.5),  which  holds  for  a  two-state  atom  (but  not  for  a 

g 

degenerate  two-level  atom).  In  Eq.  (13.5)  we  have  introduced  the  quantities 
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A®  -  A(1  +  y^)  .  (13.6) 

A®  -  Ay^  ,  (13.7) 


for  T  —  -1,  0,  1,  where 


3ire  ^ 

o 


(13.8) 


is  the  Einstein  coefficient  for  spontaneous  decay  of  this  degenerate  two-level 
atom  in  empty  space.  Furthermore,  we  write 


for  f  -  0 


for  r  -  ±  1 


(13.9) 


6  2, 

The  quantities  A  are  A*’  are  identical  in  form  to  the  relaxation  constants  A 

and  A  from  Eqs .  (10.6)  and  (10.7),  respectively,  but  their  interpretation  is 
S 

slightly  different,  as  we  shall  see  later  on. 

The  right-hand  side  of  Eq  .  (13.5)  is  probably  the  most  compact 

representation  of  the  operator  F.  As  the  first  step  towards  the  evaluation  of 

its  matrix  elements,  we  substitute  the  expressions  for  d^  and  d^ ,  and  use  the 

fact  that  a  Clebsch-Gordan  coefficient  (1  m  Irli  m  )  can  only  be  nonzero  if 

'-■g  g  '-^e  e'  ^ 

m  +  r  -  m  .  In  terms  of  the  new  constants 
g  e 


(13.10) 
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A  -  )  A®(j  m  Irlj  m  , 

“g  L  ^  g  g  '  e  e 

m  r 
e 


(13.11) 


the  operator  F  attains  the  form 


Fff  -  h  )  A  (  I  i  m  Xi  m  |<t  +  ali  m  Xj  m  I  ) 
/  m  '-^e  e  -^e  e'  '-^e  e  "^e  e' 


m 


+  W  )  A  I  I  i  m  X  j  m  1 17  +  a  I  j  m  X  j  m  I  ) 

Z  g  g  g  g  g  g  g  g 

m  ^ 
g 


(j  m  Irlj  m  )  ( j  ra'  Ir  j  j  m'  ) 
■^g  g  '-"e  e  -'g  g  '-'e  e 


m  m  r 

e  g 

m'm' 
e  g 


X  (A  <j  m  lalj  m'>lj  m  Xj  m' 1 
r-'ee'  '"^ee  gg  gg 


+  A®<j  m  IctM  ni'>|j  mXi  m'l) 
’■  g  g  g  g  e  e  -^e  e' 


(13.12) 


which  allows  a  direct  evaluation  of  matrix  elements  of  Fct  in  terms  of  matrix 


elements  of  a. 
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XIV.  EQUATION  OF  MOTION 

With  Eq.  (13.12)  we  can  expand  the  equation  of  motion  (2.6)  for  p(t)  in 
matrix  elements  This  gives  the  linear  set  of  coupled  differential  equations 


■7“  <1  m  p  i  m'>  -  -h  (A 
dt  ''ee'  ''^ee  m 


+  A  , )<j  m  pi  m'> 
m  -’e  e'  '-’e  e 
e  e 


Af(j 


m  Irh  m  )  ( j  m'lrlj  m' )<j  m  |p|i  m'> 
g  g  '-'e  e  -Jg  g  '-’e  e'  -’g  g 


m  m'  r 
g  g 


(14.1) 


T“  <i  ni 
dt  -'g  g' 


j  > 

g  g 


-h  (A 
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A  ,)<j  m 
“g  «  ® 
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m’> 
g  g 


m  m'  f 
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(14.2) 


77  <1  m  Iplj  m>-  -{h(A  +  A  )  +  iw  )<j  m  Iplj  m  >  , 

dt  ''ee'  '-^gg  m  m  o  -’ee'  '-'gg 

®  g 


(U.3) 


7~<imlp|im>--(^(A  +A  )-iw)<jmlp|im> 
dt  -'g  g  e  e  mm  o  g  g'  '-’e  e 

®  g 


(14.4) 


As  a  first  observation,  we  notice  that  Eqs .  (14.3)  and  (14.4)  for  the 


coherences  between  an  excited  state  and  a  ground  state  do  not  couple  with  the 


rest  of  the  set.  The  solution  of  Eq .  (14.3)  is  identical  to  the  solution 


(12.6)  for  a  two-state  atom,  provided  that  we  replace  A  +  A  by  A  +  A  . 

^  e  g  m  m 

&  e  g 

For  the  decay  of  the  coherence  between  two  magnetic  substates  of  the  excited 


level,  we  have  Eq .  (14.1)  with  m  m' ,  and  it  follows  that  the  time  evolution 

e  e 

of  <1  m  |p|i  m'>  is  influenced  by  the  decay  of  the  coherences  <i  m  Iplj  n>'> 
-'ee''''ee  ■’  ^  gg  gg 

between  the  various  substates  of  the  ground  level.  Note  that  the  condition 
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makes  the  term  with  m  -  m'  under  the  summation  sign  equal  to  zero, 

and  therefore  also  the  time  evolution  of  these  coherences  is  independent  of 

the  time  evolution  of  the  populations.  The  same  conclusions  hold  for  the 

coherences  in  the  ground  level,  as  described  Eq.  (14.2)  with  m  m' 

g  g 

Vhen  we  set  m'  -  m  in  Eq.  (14.1)  and  m'  -  m  in  Eq.  (14.2),  and  then 
®  ®  g  g 

perform  the  summations  over  m^  and  m' ,  resnectively ,  we  get 

^<j  m  l^lj  m>-  -A  <1  m  Iplj  m> 
dt  ■'ee'  '-'ee  m  -^ee'  '-^ee 
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m  > 
g  g 
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g  g 


IpI. 


m  > 
g  g 
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m  r 
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m  > 
e 


(14.5) 


(14.6) 


This  set  only  contains  the  populations  of  the  various  substates,  and  it  has 

the  form  of  a  set  of  rate  equations.  The  first  term  on  the  right-hand  side  of 

Eq .  (14.5)  accounts  for  the  loss  of  population  of  |j^m^>,  due  to  transitions 

to  the  various  ground  states,  and  every  m  -term  in  the  summation  gives  the 

g 

gain  of  population  of  |j^m^>,  due  to  stimulated  excitations  to  this  level  and 

from  this  particular  ground  state.  This  is  in  complete  analogy  to  Eq.  (10.10) 

for  a  two-state  atom.  The  same  interpretation  holds  for  Eq.  (14.6),  but  with 

the  role  of  | j  m  >  and  |j  m  >  reversed.  Now  that  we  have  identified  the 

®  e  g  g 

structure  of  Eqs .  (14.5)  and  (14.6),  we  can  interpret  the  significance  of  the 

relaxation  parameters  A  ,  A  ,  A®  and  A®.  From  Eq .  (14.6)  it  follows  that 

e  g 
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V  e  2 

)  A  (j  m  Ir  j  m  )  is  the  rate  constant  for  the  transition  li  m  >  -►  li  m  >. 
/^r-'gg'-’ee  ee'-^gg’ 

T 

and  the  transition  rate  equals  this  constant  times  the  population  of  |j  m  >. 
The  appearance  of  the  Clebsch-Gordan  coefficient  signifies  the  dipole 
selection  rule.  Uhen  we  sum  this  rate  constant  over  m  we  get  A  ,  as  defined 

g 

by  Eq.  (13.10),  and  from  Eq .  (14.5)  we  see  that  A  must  equal  the  rate 

m 

e 

constant  for  the  loss  of  population  of  sit'd  independent  of  which  ground 

state  this  population  ends  up  in.  This  makes  the  interpretation  of  both  rate 
constants  consistent.  Similarly,  )  A®(j  m  Irlj  m  )^  is  the  rate  constant  for 

Z  ^  &  ''’e  e^ 

T 

the  li  m  >  -•  li  m  >  transition,  and  summed  over  all  values  m  of  the  final 
g  g  e  e  e 

state  this  gives  the  rate  constant  A  for  the  loss  of  population  of  li  m  >. 

“g  S  g 

It  is  interesting  to  sum  Eqs  .  (14.5)  and  (14.6)  over  m  and  m  , 

e  g 

respectively.  This  yields 
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where 


n  -  )  <i  m 
a  l_  -'a  ( 


O  '  Q  Of 


a  -  e,g 


m 


(14.9) 


are  the  populations  of  the  two  levels,  independent  of  how  this  population  is 
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distributed  over  the  various  substates.  Only  the  relaxation  constants  A  and 

m 

e 

A  appear  in  Eqs .  (14.7)  and  (14.8),  as  is  consistent  with  our  earlier 
“g 

interpretation  of  these  parameters.  It  also  follows  immediately  that  n  (t)  + 
n  (t)  is  independent  of  time,  as  it  should  be. 


XV.  MISCELLANEA 

A.  Sum  rules  and  symmecry 

With  the  well-known  properties  of  Clebsch-Gordan  coefficients,  we  derive 
from  Eq.  (13.10) 


2j. 


1  ,e  2 
,  Aj_  +  A„ 


(15.1) 


where  the  and  f  notation  is  related  to  r-values  as  in  Eq.  (13.9).  Since 

there  are  2j^  +  1  values  of  m^  and  three  values  of  r,  the  first  equality 

states  that  the  average  value  of  A  equals  the  average  value  of  A  .  Then  the 

e 

second  equality  expresses  that  there  are  twice  as  much  parallel  degrees  of 

freedom  (x  and  y)  for  the  orientation  of  the  dipole  moment  than  there  are  for 

the  perpendicular  direction  (only  z) .  When  we  sum  A  over  m  we  find 

m  g 


irVrl  %  ''f  ■  <15-2) 

®  m  ®  r 

g 


Here , 

shall 

(15.2) 


the  left-hand  side  is  not  the  average  value  of  A  when  1  ^  .  We 

°  m  e  g 

s 

show  in  Sec.  XV. C  why  the  factor  in  front  of  the  summation  sign  in  Eq. 

must  be  l/(2j  +  1),  and  not  l/(2j  +  1). 

®  g 
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Another  intersting  relation,  which  follows  directly  from  the  properties 
of  Clebsch-Gordan  coefficients,  in  combination  with  Eqs .  (13.10)  and  (13.11), 
is 


A  -  A  ,  A  -  A  (15 . 3) 

-m  m  -m  m 

e  e  g  g 

It  can  be  shown  [31]  ,  that  the  sum  rules  (15.1)  and  (15.2)  and  the  symmetry 
relations  (15.3)  are  entirely  geometrical.  They  can  be  derived  from  the 
symmetry  of  the  vacuum  field  for  rotations  about  the  z-axis  and  reflections  in 
a  plane  through  that  axis.  Since  these  symmetries  hold  quite  generally  near  a 
flat  surface,  these  results  are  independent  of  the  details  of  the  four-wave 
mixing  process. 


B.  Fluorescence 

For  a  degenerate  two- level  atom,  the  number  of  emitted  photons  per  unit 
of  time  in  the  z  >  0  direction  equals  [19] 


1  ^ 

dt 
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h  A  Trp(t) 


)  (d  d^  +  2y  d^d  ) 
i  T  T  ■'r  T  r 


(15.4) 


When  we  substitute  the  expressions  for  d  and  d^ ,  evaluate  the  trace,  and  use 

r  T 

the  definition  of  A  ,  then  Eq.  (15.4)  reduces  to 
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(15.5) 
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where  is  Che  total  population  of  the  excited  level.  We  find  again  that  the 

first  contribution,  hAn^ ,  is  ordinary  fluorescence,  which  is  emitted 

independent  of  the  presence  of  the  PC.  Every  term  under  the  summation  sign 

corresponds  to  the  occurrence  of  a  three-photon  process  between  the  particular 

state  (j  m  >  and  the  excited  level.  It  appears  that  the  contributions  from 
S  & 

the  various  excited  states  combine  in  the  single  rate  constant  A  ,  and 

m 

g 

therefore  the  rate  of  emission  of  fluorescent  photons  due  to  a  population  of 

li  m  >  is  equal  to  the  rate  of  spontaneous  transitions  from  |i  m  >  to  the 
'  g  g  '  g  g 

excited  level,  as  could  be  expected. 


C.  Ideal  PC 


For  the  case  of  a  degenerate  two-level  atom  and  a  transparent  PC,  we  can 


find  the  limit  of  an  ideal  PC  by  setting 


yi-.y|j-mpl2 


g  a 

With  Eqs.  (13.6)  and  (13.7)  we  then  obtain  A^  and  A^,  which  become  independent 
of  r .  In  Eqs.  (13.10)  and  (13.11)  we  can  then  perform  the  summations,  and  it 


follows  that  A  and  A  become  independent  of  m  and  m  ,  respectively.  Hence 
e  “g  eg 

we  define  A  and  A  as 
e  g 


A 

e 


(15.7) 


and  they  are  found  to  be 
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A  -  A®  -  A(1  +  ^ 
e  T 


(15.8) 
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Equation  (13.5)  reduces  to 
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(15.10) 


in  terms  of  the  projectors 


P  -  )  |j  m  Xj  ml  ,  a  -  e,g 
m 

a 


(15.11) 


For  the  populations  we  now  find 

(15.12) 

(15.13) 

which  is  identical  in  form  to  Eqs.  (10.10)  and  (10.11).  This  shows  that 
and  A^  are  the  rate  constants  for  the  depopulation  of  the  excited  level  and 
the  ground  level,  respectively,  Irrespective  of  how  this  population  is 
distributed  over  the  various  states.  Notice  that  the  excited  level  and  the 
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ground  level  have  the  same  relaxation  constant  as  their  individual  states 

|j  m  >  and  li  m  >,  as  follows  from  Eq.  (15.6).  It  can  also  be  shown  that  the 
Ue  e  g  MV/ 

coherences  between  the  two  levels  decay  with  h(A  +  A  )  as  inverse  relaxation 

eg 

constant,  and  that  the  fluorescence  rate  educes  to 


Z - H  An  +  A  n 

dt  egg 


(15.14) 


We  conclude  that  in  the  limit  of  an  ideal  PC  this  degenerate  two- level  atom 

behaves  similarly  as  a  two-state  atom,  at  least  as  far  as  the  time  evolution 

of  the  atomic  density  operator  is  concerned,  or  the  fluorescent  yield.  Notice 

that  also  the  expression  (10.6)  and  (15.8)  for  the  relaxation  constant  of  the 

excited  state  are  identical  in  form,  but  that  A  in  Eq.  (15.9)  acquires  an 

g 

additional  factor  (2j  +  l)/(2j  +1),  as  compared  to  A  in  Eq.  (10.7).  This 

®  g  g 

factor  equals  the  ratio  of  the  number  of  excited  states  and  the  number  of 
ground  states,  and  it  also  appears  explicitly  in  expression  (15.10)  for  F  and 
in  the  sum  rule  (15.2)  for  the  average  of  A  It  takes  into  account  the 

IQ 

g 

difference  between  the  number  of  excited  states  and  the  number  of  ground 
states,  and  it  leads  to  the  correct  form  of  the  set  of  rate  equations 
(15. 12), (15. 13). 


XVI.  CONCLUSIONS 

We  have  studied  the  relaxation  of  an  atom  in  the  vicinity  of  a  PC.  It 
appeared  that  the  atomic  relaxation  operator  F  could  be  expressed  in  the  two 
vacuum  correlation  functions  fj^(w)  and  f||(a>),  which  contain  all  necessary 
information  about  the  medium  and  the  four-wave  mixing  process.  We  were  able 
to  express  these  functions  in  terms  of  the  Fresnel  reflection  and  transmission 
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coefficients  for  plane  waves,  and  the  normal  distance  between  the  atom  and  the 

surface  of  the  medium.  These  coefficients  depend  in  a  complicated  way  on  both 

the  linear  and  nonlinear  interactions,  the  polarization  and  the  angle  of 

incidence.  For  instance,  it  was  shown  in  a  recent  experiment  [39]  that  even 

Che  Fresnel  coefficient  for  the  reflection  of  the  specular  wave  is  strongly 

affected  by  the  nonlinear  interaction.  We  found  that  the  correlation 

functions  could  be  expressed  in  terms  of  the  parameter  functions  x  ,  y  and 

a  ■'  a 

z  ,  with  a  -  or  ||  . 
a  " 

Subsequently  we  worked  out  the  model  case  of  a  two-state  atom,  in 
combination  with  the  approximation  of  an  ideal  PC.  The  relaxation  operator  F 
then  assumes  the  simple  form  (10.5).  We  obtained  the  relaxation  constants  A 

e 

and  A  for  the  excited  state  and  the  ground  state,  respectively,  and  we 
S 

identificed  the  four  processes  which  are  responsible  for  the  relaxation  (Fig. 

2)  .  Ordinary  spontaneous  fluorescent  emission  in  combination  with  an  |e>  -♦ 

|g>  transition  gives  the  usual  contribution  to  A^  (equal  to  A).  Furthermore, 

stimulated  |e>  -*  | g>  transitions  in  combination  with  the  emission  of  a  photon, 

and  stimulated  three -photon  processes  which  excite  the  atom  contribute  to  the 

relaxation.  The  rate  constants  for  both  (down  and  up)  stimulated  processes 

are  identical,  at  least  in  the  limit  of  an  ideal  PC  in  combination  with  a  two- 

state  atom.  This  could  be  expected  from  the  analogy  with  stimulated  thermal 

transitions,  although  the  physical  mechanism  is  different.  These  stimulated 

processes  always  increase  the  values  of  the  relaxation  constants,  thereby 

decreasing  the  relaxation  time  of  the  system.  The  three-photon  processes  are 

responsible  for  a  finite  population  of  the  excited  state  in  the  limit  t  ®. 

We  found  that  both  A  and  A  are  independent  of  the  atom- surface  distance, 

eg 

which  resembles  similar  predictions  by  Agarwal  [40]  and  Cook  and  Hilonni 


[41]. 
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FIGURE  CAPTIONS 

Fig.  1  ■  A  plane  wave  is  incident  upon  a  layer  of  nonlinear  material  with 
thickness  A.  At  the  boundary  z  -  0  this  wave  partially  reflects  into  the 
specular  direction  (r-wave),  and  partially  travels  into  the  medium.  At  z  -  -A 
there  is  again  partial  reflection  and  partial  transmission,  which  gives  rise 
to  multiple  reflections  irside  the  medium  and  the  t-wave  in  z  <  -A.  These 
waves  are  indicated  by  the  solid  arrows,  where  the  arrow  head  indicates  the 
propagation  direction.  These  waves  are  already  present  for  a  linear  medium. 
Due  to  the  nonlinear  interaction,  the  two  waves  inside  the  medium  couple  each 
with  a  counterpropagating  negative -frequency  wave  (provided  that  the  incident 
field  has  a  positive  frequency),  which  leave  the  medium  at  the  boundaries  z  - 
0  and  z  -  -A  as  the  pc -wave  and  the  ni-wave,  respectively.  These  negative- 
frequency  waves  are  shown  as  broken  arrows.  For  a  dielectric  constant  equal 
to  unity,  the  r-wave  and  the  ni-wave,  including  the  two  corresponding  waves  in 
the  medium,  disappear  to  a  good  approximation.  The  figure  illustrates  the 
situation  for  which  the  incident  wave  comes  from  the  positive  z-direction.  As 
explained  in  the  text,  there  are  also  vacuum  waves  incident  from  the  region 
z  <  -A,  and  they  give  rise  to  ni-waves  and  t-waves  in  z  >  0. 

Fig.  2.  Illustration  of  the  four  processes  which  contribute  to  the  atomic 
relaxation  and  to  the  emission  of  observable  fluorescence.  The  rates  for  the 
processes  are  indicated,  and  the  wlggly  arrows  show  pictorially  the 
propagation  direction  of  the  photons  involved.  Diagram  (1)  represents 
ordinary  spontaneous  decay,  together  with  an  |e>  -►  |g>  transition  of  the  atom. 
The  fluorescent  photon  travels  into  the  direction  of  the  detector  and 
contributes  to  the  fluorescent  yield,  whereas  the  atomic  decay  gives  rise  to 
relaxation.  The  same  process  is  shown  in  diagram  (2),  except  that  the  photon 
here  travels  towards  the  PC.  Therefore,  this  process  contributes  to  the 
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relaxation,  but  not  to  the  observable  fluorescence.  Dlagraa  (3)  looks  the 
same  as  diagram  (1),  except  that  the  rate  is  multiplied  by  |p]  We  identify 
this  process  as  stimulated  decay  combined  with  the  emission  of  a  photon.  This 
photon  also  travels  into  the  positive  z-direction,  because  it  is  a  t-wave,  but 
it  does  not  contribute  to  the  fluorescence.  In  the  three -photon  process  from 
diagram  (4)  the  atom  makes  the  transitions  |g>  -►  |  e>  |  g>  where  the 
first  and  the  last  transitions  are  accompanied  by  the  absorption  of  a  photon 
with  frequency  u.  During  the  decay  in  between  these  absorptions,  the  atom 
emits  a  photon  with  frequency  in  the  positive  z-direction.  Notice  that 
emitted  photons  in  a  stimulated  transition  always  travel  in  the  positive  z- 
direction. 


